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PREFACE 


In  January  1985  the  Director  of  the  Science  and  Technology  Directorate  of  the 
Strategic  Defense  Initiative  Organization  (SDIO)  asked  members  of  the  IDA  research  staff 
to  investigate  the  feasibility  of  developing  a  y-ray  laser.  The  staff  determined  what  work 
had  been  done,  who  was  currently  working  in  the  field,  and  what  work  should  be 
encouraged  or  supported.  A  workshop  was  convened  for  researchers  directly  involved 
both  in  gamma-ray  laser  work  and  in  ancillary  fields  such  as  nuclear  structure,  radiation 
propagation  in  crystals,  Mdssbauer  Effect,  and  q>tical  lasers.  Next,  an  in-house  study  was 
undertaken  to  clarify  critical  issues  concerning  the  various  pumping  schemes  proposed  at 
the  workshop  as  well  as  systems  questions  about  the  y-ray  laser  as  a  working  device. 

The  proceedings  of  the  workshop  were  published  in  the  form  of  a  report  to  the 
Irmovative  Science  and  Technology  Office  (1ST)  of  the  SDIO.  The  work  completed  in 
1985  was  presented  in  IDA  Paper  P-2021. 

In  1986,  the  in-house  work  focused  on  extending  the  data  base,  the  nature  of 
superradiance  in  the  y-ray  laser  context,  and  a  detailed  investigation  of  the  upconversion 
pumping  scheme.  A  discussion  of  nuclear  systematics,  investigations  of  electron-nuclear- 
driven  pumping,  and  lifetime  measurements  rounded  out  that  study.  The  results  of  the 
FY1986  effort  were  presented  in  IDA  Paper  P-2004. 

In  1987,  the  IDA  research  staff  looked  at  the  state  of  the  art  and  assessed  the 
situation  in  y-ray  laser  work  with  focus  on  two  areas  of  research  interest  critical  to  concepts 
for  developing  a  y-ray  laser.  Heating  effects  associated  with  upconversion  techniques  were 
discussed.  The  sources  of  inhomogeneous  broadening  which  destroy  the  Mdssbauer 
Effect  were  investigated  and  techniques  available  for  restoring  the  resonance  by  external  or 
internal  fields  were  ccmsidered  (IDA  Paper  P-2083). 

I  In  1988,  the  in-house  work  concentrated  on  (1)  establishing  a  theoretical  model  of 

nuclear  superfluorescence  that  took  into  account  specific  characteristics  of  nuclear  radiation 
emission  and  transport  that  are  important  for  superfluorescence  on  the  nuclear  scale, 
(2)  examining  in  more  detail  and  for  more  realistic  parameters  the  heating  effects  inherent  in 


iii 


upconversion  concepts,  and  (3)  examining  the  speed  of  response  of  nuclei  to  applied 
external  fields,  as  exhibited  by  their  spectra  (IDA  Paper  P-2175). 

In  1989  the  work  concentrated  on  developing  a  OKxlel  of  nuclear  superfluorescence 
that  included: 

1 .  Multilevel  systems, 

2 .  Incoherent  pumping  at  finite  rates, 

3.  Spatial  attenuatiai  of  the  electromagnetic  field  in  the  active  region, 

4 .  Dephasing  due  to  inhomogeneous  teoadening, 

5 .  Population  depletion  due  to  competing  processes  like  internal  conversion  and 
spontaneous  transverse  emission. 

This  was  done  by  generalizing  the  Haake-Reibold  model  of  superfluorescence  based 
on  Maxwell-Bloch  equations.  The  effects  of  inhomogeneous  broadening  on 
superfluorescence  were  studied  using  a  theory  developed  by  J.  Ebeiiy.  The  enhancement 
of  superfluorescence  in  inhomogeneously  broadened  systems  time-dependent  hyperfine 
interactions  was  investigated  using  the  theory  developed  by  M.  Blume.  We  also 
investigated  coherent  aspects  of  nuclear  emission  and  absorption  in  the  presence  of 
homogeneous  and  inhomogeneous  broadening  of  source  and  absorber  (IDA  Paper 
P-2335). 

The  present  document  describes  the  1990  work  in  which  we  iq}plied  the  Maxwell- 
Bloch  model  (and  other  models)  of  the  interaction  of  the  electromagnetic  field  with  matter  to 
the  y-ray  laser  problem  in  the  following  way: 

1 .  We  compared  amplified  spontaneous  emission  (ASE)  and  superfluorescence 
(SF)  within  the  same  model  and  found  the  conditions  for  SF  to  be  more  easily 
attainable. 

2.  We  analyzed  the  feasibility  of  SF  in  a  real  nucleus,  ^^o  assumed  prepared  in 
the  isomeric  (inverted)  state  through  a  thermal  neutron  reaction  with  the  parent, 
and  calculated  SF  pulse  intensity  under  different  conditions.  This  is  the  first 
published  result  of  a  real  nucleus,  pumped  through  an  existing  mechanism  and 
achievable  fluxes  of  the  pumping  radiation,  to  give  a  SF  pulse. 

3.  From  the  Monte  Carlo  runs  with  the  Maxwell-Bloch  equations  we  developed 
parametric  equations  characterizing  the  SF  emission  delay  time.  The 
parametric  equations  simplify  future  analyses  of  SF  phenomena. 

4.  We  tested  the  popular  assumption  that  SF  in  an  attenuating  medium  can  be 
modeled  by  assuming  a  reductitm  in  the  effective  length  of  the  cylindrical  active 


volume,  in  effect  assuming  the  number  of  collective  radiators  to  be 
proportional  to  the  inverse  of  the  attenuation  coefficient.  This  assumption 
proved  to  be  inaccurate  and  unacceptable  for  candidate  selection. 

The  development  of  a  y~ny  laser  is  viewed  as  a  high-risk/high-payoff  undertaking. 
IDA'S  involvement  focuses  on  minimizing  the  risk  and  on  striving  to  redirect  the  effort 
when  proposed  schemes  are  shown  not  to  be  feasible. 
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ABSTRACT 


This  paper  compares  two  processes  that  can  occur,  under  different  conditions,  in 
the  same  system  of  nuclei:  amplified  spontaneous  emission  (ASE)  and  superfluoiescence 
(SF).  It  concludes  that  the  (Trammell  and  Hannon)  conditions  for  the  occurrence  of  ASE 
are  more  restrictive  than  those  for  the  occurrence  of  SF.  Therefore,  the  conclusion  by 
some  authors  based  on  an  analysis  of  ASE,  that  cohoent  electromagnetic  emission  from  a 
nuclear  system  is  not  feasible,  is  incorrect  The  paper  examines  the  feasibility  of  nuclear 
SF  in  It  concludes  that  SF  will  occur  in  the  excited  state  obtained  by  thermal  neutron 
pumping  with  fluxes  ranging  from  10^^  to  10^  neutrons  per  second. 
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SUMMARY 


This  report  describes  the  1990  research  effort  by  members  of  the  IDA  staff  in  the 
field  of  y-ray  lasers.  The  work  is  part  of  a  continuing  task  in  support  of  the  Innovative 
Science  and  Technology  Office  (1ST)  of  the  Strategic  Defense  Initiative  Organization 
(SDIO).  The  development  of  a  y-ray  laser  is  a  high-risk  science  and  technology 
undertaking.  IDA  involvement  has  focused  in  large  measure  on  minimizing  the  risk  and 
attempting  to  redirect  the  program  as  quickly  as  possible  when  proposed  schemes  prove 
infeasible. 

In  this  report  we  look  for  optimal  and  realistic  conditions  for  the  emission  of 
coherent  electromagnetic  radiation  from  nuclei  and  the  development  of  a  y-ray  laser.  We 
compare  the  requirements  for  amplified  spontaneous  emission  (ASE),  or  pulsed  laser 
emission,  with  superfluorescence  (SF),  a  cooperative  emission  from  an  inverted  population 
of  nuclei.  We  find  the  requirements  for  SF  less  restrictive  on  the  nuclear  and  solid  state 
properties  of  the  acdve  medium.  Finally  we  use  the  general  theory  of  SF  to  examine  the 
feasibility  of  nuclear  SF  in  obtained  in  the  excited  state  by  thermal  neutron  pumping 
with  fluxes  ranging  from  10^*^  to  1023  neutrons/sec.  It  should  be  emphasized  that  our 
calculation  is  based  on  real  properties  of  materials  and  that  this  is  the  first  calculation 
published,  to  our  knowledge,  that  includes  broadening  effects,  attenuation,  coupling 
transitions  and  finite  pumping  times  in  a  single  model. 

The  interaction  of  the  radiation  field  with  matter  leads  to  coherent  emission  (ASE) 
and  cooperative  emission  (SF)  in  the  same  inverted  amplification  medium.  The  two 
phenomena  can  be  distinguished  theoretically  and  experimentally.  The  selection  of 
appropriate  candidate  nuclei  depends  on  the  specific  phenomena  because  of  the  different 
requirements  that  have  to  be  met.  We  found  that  the  SF  process  for  nuclei  is  less 
restrictive,  whereas  researchers  have  generally  used  the  Schawlow-Townes  threshold 
condition  for  ASE  as  a  guideline  for  selecting  candidates.  The  Schawlow-Townes 
condition  is  derived  in  the  linear  approximation  whereas  we  use  the  full  nonlinear 
interaction  Hamiltonian  to  obtain  our  results. 

We  concentrated  on  the  SF  process  and  studied  the  requirements  for  SF  pulse 
emission  as  a  function  of  the  number  of  cooperating  nuclei,  the  attenuation  properties  of  the 

S-1 


medium,  the  natural  single  nuclear  lifetime  and  the  pumping  rate.  We  characterize  the 
emitted  pulses  by  the  delay  time  and  the  intensity  and  use  the  results  from  the  Monte  Carlo 
solutitms  of  the  Maxwell-Bloch  equations  to  obtain  parametric  equations  for  the  delay  time. 
This  allows  us  to  make  predictions  and  perform  certain  analyses  much  more  easily  without 
additional  time-consuming  Monte  Carlo  calculations. 

A  popular  approach  for  taking  into  account  the  attenuation  of  the  medium  has  been 
to  assume  a  leductitm  in  the  active  volume  by  the  inverse  of  the  attenuation  coefficient  We 
checked  this  assumption  by  performing  the  complete  Monte  Carlo  calculation  of  the 
Maxwell-Bloch  equations  and  found  a  discrepancy  of  a  factor  of  four  in  the  delay  time 
calculations.  The  Monte  Carlo  calculation  gave  shorter  delay  times  than  the  volume 
reduction  model.  This  difference  could  lead  to  a  critical  error  in  the  candidate  selection 
process  and  rule  out  certain  nuclei  that  could  in  fact  be  useful. 
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I.  INTRODUCTION 


In  the  history  of  y-ray  laser  research,  initially,  attention  focused  on  the  search  for 
concepts  and  nuclei  that  satisfied  the  Schawlow-Townes  (S-T)  condition  for  lasing,  which 
can  occur  either  in  the  form  of  continuous  wave  (CW)  or  pulsed  emission.  Lasing  in  the 
form  of  pulse  emission  initiated  by  spontaneous  decay  processes  is  usually  referred  to  as 
amplified  spontaneous  emission  (ASE).  Recognition  that  in  the  nuclear  regime  ASE  is  a 
more  realistic  phenomenon  to  expect  than  CW  lasing  changed  the  direction  of  the  Y*ray 
laser  research. 

Trammell  and  Hannon  (Ref.  1)  analyzed  the  conditions  required  for  the  occurrence 
of  ASE  and  found  that  it  is  a  more  restricted  process  than  lasing  and  is  limited  by  a 
modified  version  of  cooperative  emission  following  inversion,  the  S>T  condition.  They 
also  drew  a  similar  conclusion  about  the  phenomenon  of  superfluorescence  (SF);  however, 
reasons  for  disagreeing  with  that  inference  are  discussed  in  Chapter  m. 

More  recently  SF  was  considered  a  more  feasible  and  promising  approach  to 
nuclear  lasing  (Ref.  2).  In  subsequent  investigations  the  processes  of  ASE  and  SF  have 
been  treated  in  more  detail  and  generality.  The  transition  region  between  ASE  and  SF  has 
also  been  investigated  experimentally  in  atomic  systems  (Ref.  3a,  b).  Such  studies  provide 
the  background  and  permit  compariscms  of  ASE  and  SF  in  nuclear  systems  leading  to  valid 
checks  on  candidates  for  experimentally  demonstrable  feasibility  of  ASE  or  SF. 

Nevertheless,  not  all  of  the  important  physical  effects  that  should  be  taken  into 
account  in  the  analysis  of  real  nuclear  systems  have  been  considered  in  a  single  theory. 
The  most  complete  theory  for  treating  nuclear  SF  was  presented  in  IDA  paper  P-2335  and 
is  used  in  our  analysis. 

In  Chapter  n  we  compare  the  necessary  conditions  for  the  existence  of  the 
phenomena  ASE  as  modeled  by  Trammell  and  Hannon  and  SF  and  show  that  the 
occurrence  of  SF  is  more  probable  for  the  nuclear  case  under  certain  conditions. 

Chapter  m  deals  with  the  details  of  a  parametric  description  of  the  SF  processes, 
(derived  from  the  Monte  Carlo  calculations  given  in  Ref.  4).  This  approach  permits 
investigations  of  SF  phenomena  without  resort  to  the  lengthy  calculations  of  Ref.  4. 
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Chapter  IV  uses  the  earlier  theoretical  development  to  investigate  the  feasibility  of 
SF  in  a  real  nucleus,  ^Co. 
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11.  COMPARISON  OF  AMPLIFIED  SPONTANEOUS 
EMISSION  AS  MODELED  BY  TRAMMELL  AND 
I  HANNON  AND  SUPERFLUORESCENCE 


I 


A.  LASING  (AMPLIFIED  SPONTANEOUS  EMISSION) 


ASE  can  be  viewed  as  a  transptnt  process  with  gain.  It  is  governed,  as  is  lasing  in 
general,  by  the  Schawlow-Townes  threshold  condition  (Ref.  1)  which  simply  states  that 
for  gain  to  occur  in  a  region,  the  stimulation  of  photons  has  to  exceed  the  absorption.  For 
this  to  occur,  the  gain  coefficient  has  to  be  positive; 

K  =  OsAn*-p.>0  , 


where  the  stimulation  cross-section  0$  is  given  by 


”  2n  (1  +  a)  (1  +  a) 


(1) 


and 


» 


I 


I 


An*  =  inversion, 

p  =  linear  attenuation  coefficient, 

P  =  branching  ratio, 

X  s  the  wavelength  of  the  emission, 

a  =  inhomogeneous  broadening  parameter, 

a  =  internal  conversicm  coefficient, 

f  =  recoilless  fraction. 

The  Schawlow-Townes  condition  sets  the  threshold  for  the  occurrence  of  both 
lasing  and  ASE.  The  basic  process  is  a  transport  process  independent  of  tinae  (this  is 
strictly  true  only  for  a  steady  state  operation).  For  a  pulsed  laser  or  ASE  the  conditions 
must  allow  for  the  fact  that  the  population  of  the  excited  level  decreases  during  emissitm. 
Only  the  inversion  An*  changes;  the  other  parameters  in  the  Schawlow-Townes  condition 
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are  independent  of  time.  This  is  discussed  by  Trammell  and  Hannon  in  Ref.  1.  where  the 
threshold  condition  for  pulsed  lasers  is  derived.  We  use  their  results  in  our  comparison  of 
ASE  and  superfluorescence  (SF).  We  also  extend  their  calculations  into  regions  of 
particular  interest  and  present  the  results  in  Appendix  A. 

B.  SUPERFLUORESCENCE 

SF  can  be  viewed  as  a  spontaneous  emission  by  coupled,  cooperating  radiators. 
The  process  is  initiated  by  quantum  fluctuations  of  the  vacuum  electromagnetic  field,  or  the 
emission  of  a  photon,  which  gives  rise  to  a  phasing  or  build  up  of  correlations  in  the 
inverted  region.  SF  is  governed  by  the  dynaimcs  of  the  evolving  system  as  described,  for 
example,  by  the  Maxwell-Bloch  equations  with  random  noise  soinoes  (Ref.  4,  Chapter  I). 
No  simple  threshold  condition  connecting  the  pertinent  parameters  exists  as  in  the  case  of 
lasing.  Conditions  for  the  occurrence  of  SF  are  described  in  terms  of  characteristic  times 
representing  the  speed  of  competing  processes:  namely,  the  single  nucleus  lifetime  Tn*  the 
radiative  lifetime  Xr  the  cooperation  time  Tc,  the  time  of  eirussion  of  the  SF  pulse  XsF,  the 
delay  time  of  the  pulse  xd>  ttnd  the  dephasing  time^  and  the  pumping  time  ty.  The 
cooperation  time^  Xq  sets  the  absolute  upper  limit  on  the  SF  time  tsF  and  consequently  on 
the  delay  time  td-  We  discuss  this  in  the  next  section.  Intuitively  one  expects  that  for  SF 
to  occur  the  following  conditions  should  hold: 

“ty  <  Xjj  <  X^,  X||,  Xp  . 

More  rigorous  analysis  (Ref.  5)  shows  that 


is  actually  the  appropriate  condition  between  the  delay  time,  SF  time,  and  the  dephasing 
time. 

^  This  commonly  is  due  to  collisions  or  inhomogenous  tnoadening. 

2  The  cooperation  time  is  obtained  from  the  Dicke  formula  for  the  emission  rate,  it  =  -L  Nc  —  X2  /A, 
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whoe  Nc  is  the  number  of  emitters,  \  the  linewidih  of  the  radiation  and  A  the  cross  section  of  the  tod 
forming  the  inverted  region  and  a  self  consistent  argument  determining  how  many  atoms  can  be  covered 

by  the  emission  of  a  single  atom  during  the  Dicke  emission  time,  With  Nc=  Ap  we  get 


These  times  are  defined  in  terms  of  elementary  parameters  or  derived  from  basic 
theories  (Refs.  5, 6, 7).  For  our  purposes  they  are  as  follows: 

tr  =  ax„  . 

yjcprk 
^  X  pi 

to  =  ix5p/n(N)  . 


where  y  is  the  single  nucleus  spontaneous  pumping  rate  producing  the  inversion,  F  is  the 
single  nucleus  spontaneous  radiative  rate,  p  is  the  nuclear  density,  a  is  the  inhomogeneous 
broadening  parameter  (Ref.  1)  defined  by  Fp  ^  (1+a)  Fq^  and  /  the  length  of  the  active 
region.  The  expressions  of  tSF  and  td  given  in  Eqs.  (6)  and  (7)  are  model  dependent  and 
are  used  here  for  convenience  (see,  for  example.  Ref.  8  or  Cluq)ter  m).  The  phenomena  of 
SF  is  discussed  in  greater  detail  in  Refs.  9  and  10. 

Using  the  Maxwell-Bloch  equations  with  stochastic  sources  (Ref.  4,  Chapter  I),  we 
have  studied  the  feasibility  of  nuclear  SF.  This  formalism  characterizes  SF  under 
physically  meaningful  conditions.  Processes  important  to  nuclear  SF  such  as  attenuation, 
inhomogeneous  line  broadening  and  competing  transitions  are  all  easily  and  naturally 
incorporated.  The  formalism  also  allows  consideration  of  finite  pumping  times  and 
incoherent  pumping  mechanisms.  Using  the  delay  time,  ti>,  as  the  criterion  for  the  strength 
of  SF  emission,  we  found  the  effective  length  left  of  the  active  medium  for  dissipative 
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systems  (p  ^  0)  to  be  about  /cff = |[  (or  higher)  instead  of  -jj  as  is  often  assumed  (Refs.  1, 

11).  We  used  the  Maxwell-Bloch  equations  and  the  Eberly  theory  (Ref.  7)  to  study  the 
effect  of  inhonoogeneous  broadening,  another  phenomenon  critical  in  the  generation  of 
nuclear  SF.  Inhomogeneous  Inoadening  produces  dephasing  with  a  characteristic  time  of 
so  that  at  short  times  after  inversion  (t«T^)  the  effect  on  SF  of  inhomogeneous 
broadening  is  small  (increasing  with  time),  but  not  constant  In  the  equation  for  Os  and  the 


gain  coefficient  (the  Schawlow-Townes  condition),  the  effect  of  inhomogeneous 

broadening  is  a  constant  reduction  factor  •  'fhus  SF  and  lasing  are  quite  different  in 

this  respect;  both  the  attenuation  and  inhomogeneous  broadening  play  a  different,  and  as 
will  be  seen,  less  restrictive  role  in  the  SF  process  than  in  the  lasing  process. 

C.  LIMIT  ON  THE  COOPERATION  LENGTH 

The  absolute  upper  limit  on  effective  length  of  the  active  medium  can  be  found  from 
the  cooperation  dme.  The  cooperation  time  is  determined  by  the  number  of  nuclei  that  can 
interact  widiin  the  lifetime  of  the  nuclear  state  (Ref.  6).  From  Eq.  (2)  we  have 


where  r  =  7~  is  the  radiadve  rate. 

The  maximum  effective  length  has  to  be  less  than  the  conelation  length,  i.e., 

V 

From  this  we  can  get  the  absolute  limit  on  the  SF  dme; 


Using  typical  nuclear  parameters, 

p  =  1022  , 

X,  =  lA  =  lO^cm  , 
a  =  10  , 

c  =  3  X  10*Ocm/sec  , 
cpX,2  =  3  X  10*^  sec~i  , 


(4) 


(5) 
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and  assuming  f  =  1,  we  get 


► 


» 


I 


V^(f)^A=2.3xlO-‘/;  .  (6) 

This  gives  the  minimum  SF  pulse  emission  time.  Generally  the  experimentally 
observed  pulse  emissitm  time  is  greater  than  tsF> 


D.  COMPARISON  OF  CONDITIONS  ON  ASE  AND  SF 

For  the  steady  state  laser  the  Schawlow-Townes  gain  condition  is 

^  “  27c  (l4a)  (1+a)  • 


(7) 


As  shown  by  Trammell  and  Hannon  (Ref.  1),  this  has  to  be  modified  for  die  pulsed 
laser,  or  ASE,  to  account  for  the  reductions  in  the  inversion  during  emission.^  Their 
results  show  that  a  factor  of  21  multiplying  ^  in  the  condition  on  K  is  required,^  thus  for 
the  present  discussion  we  assume  that 

f  An*  ^ 

'^~2jt  (l+a)(l+a)~21li>0  (8) 

for  a  pulsed  laser. 

Assuming  oc»l,  a»l. 


we  get 


and  finally 


fAn*  ^  j 
2ic  aa(21)[i 


0241  < 


X^fAn* 

2ii(21)  ‘ 


3  Thuisthesaineasthewell-known"laserlethargy''ofM.  Scully  and  F.  Hoffman. 

*  This  condition  has  limited  qiplicability  to  sy^ms  where  the  decrease  in  invetsion  is  assumed  to  be 
^  due  to  the  natunl  decay  proc^.  For  Mgh-gain  systems  the  factor  of  21  increases,  as  we  discovered  in 

our  investigations  when  we  generalized  the  results  of  Trammell  and  Hannon  (see  Appendix  A). 
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The  SF  condition  requires  that  SF  dynamics  are  faster  than  dephasing  or 
depopulating  piticesses.  We  assume  •»  0  and  use 

as  the  mathematical  statement  of  this  amdititx). 

Assuming  an  effective  length  for  p  ^  0,  as  determined  by  calculations 

T 

based  on  the  Maxwell-Bloch  equations  (Ref.  4)^  using  t^=  — ,  \  =  ax„  and  applying  the 
above  criteria  we  get: 


2  tgp  In  (N)  tjpj  <  ^ 


■//nTr<^  , 


^  X^fp  (4) 


and  therefore 


<1  . 


‘2  xV 


(paa)<- 


3X^p 

jrV2/n(N) 


^  Detailed  calculations  are  given  in  Chapter  m. 


In  making  the  ctnnparison  between  SF  and  ASE  we  will  consider  the  constraints  cm 
the  product  (a  a  4)  for  the  competing  processes.  The  internal  conversion  coefficient  a  is  a 
given  for  a  particular  nuclear  isomer,  and  we  don’t  have  a  lot  of  external  control  over  the 
existing  value.^  Nevertheless,  it  is  an  important  parameter  for  candidate  selection.  The 
inhomogeneous  broadening  represented  by  a,  on  the  other  hand,  can  be  controlled  by 
external  RF  pulses  and  other  techniques,  and  the  linear  attenuation  coefficient  p.  can  be 
reduced  through  the  Borrmann  Effect  Conveniently,  the  product  (a  a  p)  appears  in  the 
requirements  for  both  ASE  and  SF  processes. 

The  constraint  (xi  ASE  is 


and  on  the  SF  process 


aap< 


3X^fp 
7c/2/n  (N) 


Taking  the  ratio  of  the  right-hand  sides  of  Eqs.  (11)  and  (12)  we  get 

p  -  / 

JcV2/n(N)  / 

which  reduces  to 

126  ^ 

V2/n(N) 

ifAn*  =  p. 

For 


N=  102, 

R  -  41.5 

N=103, 

R  »  33.9 

N=104, 

R  -  29.4 

N=  106, 

R-24.0 

N  =  10^2 

R  -  16.9 

N  =  1022 

R  »  12.5 

An* 

2Jt(21) 


89.095 

V/iiN 


(12) 


(13) 

(14) 


^  Some  reduction  can  be  achieved  by  separation  (dinner  shell  dectronfirom  the  nuclear  core. 


This  comparison  shows  that  SF  permits  larger  values  of  (a  a  p)  and  is  therefore 
less  restrictive  with  respect  to  the  nuclear  and  atomic  requirements  of  the  active  medium 
than  ASE  as  modeled  by  Tramell  and  Hannon. 

Finally,  if  we  remove  the  assumption  that  a  »  1,  a»l  and  p  >  0,  then  the 
appropriate  conditions  for  ASE  and  SF  are 

2 

p(a+l)(a  +  l)<M^  (15) 


for  ASE,  and 


(a+l)(l+a)< 


(16) 


m.  PARAMETRIC  FITS  TO  MONTE  CARLO 
SIMULATIONS  OF  SUPERFLUORESCENCE 
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We  have  studied  SF  phenomena  using  Monte  Carlo  solutions  (Ref.  4)  of  Maxwell- 
Bloch  equations.  The  details  of  the  calculations,  which  required  a  considerable  amount  of 
computer  time,  and  our  findings  are  described  in  two  previous  IDA  reports  (Refs.  4,  8). 
To  save  computer  time  in  analyzing  and  predicting  experimental  results  we  have  developed 
parametric  iiKxlels  based  on  these  Mrmte  Cario  calculations. 

Fig.  m-la  shows  the  delay  time  normalized  to  x^p  as  a  function  of  the  cooperation 
number  N.  In  the  runs  the  attenuation  coefficient  p  =  0,  the  dots  represent  Monte  Carlo 
runs  and  the  curve  represents  the  least  square  fit  to  the  parametric  equation 

XD  =  a[l +bln(N)]ln(N),  (17) 

in  which  the  parameter  values  corresponding  to  die  best  tit  are 

a  =  2.43  , 
b  =  0.02  . 

Figure  m-lb  shows  a  parametric  tit  to  the  same  data  in  terms  of  n  =  In  N,  of  the 
equation 

Xt,  =  x  T 1+P  n+pVl, 

D  DOL  2  J  (Ig) 

in  which  the  parameter  values  are: 

Too  =  0.158  , 

Pi  =67.885  , 
p2  =7.067  . 

In  their  theoretical  treatment  Bonifacio  and  Lugiato  (Ref.  9)  [also  cf.  Gross  and 
Haroche  (Ref.  10)]  obtained  a  logarithmic  dependence  on  the  cooperation  number  N,  while 
Polder,  Schuurmans,  and  Vrehen,  in  their  work  (Ref.  12)  derived  a  dependence  on  the 
square  of  the  logarithm  of  N.  Our  result,  based  on  the  work  of  Haake  and  Reibold 
(Ref.  13)  is  a  linear  combination  of  both. 
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Figure  Ill-la.  Parametric  fit  of  delay  time  to  obtained  from 
Monte  Carlo  rune  as  a  function  of  N. 


We  also  ran  Monte  Carlo  calculations  for  N  =  10^,  1(P,  10^,  and  different  values 
of  n.  The  time  delay  to  versus  p.  is  plotted  in  Fig.  111-2  for  the  different  values  of  N.  In 
Fig.  in-3  we  show  the  normalized  area  under  the  pulse  as  a  function  of  p..  Note  how 
closely  the  curves  overlap,  indicating  little  effect  of  the  cooperation  number  on  the  change 
in  total  emission  with  attenuation.  The  effect  of  cooperation  number  on  the  delay  time 
however  is  substantial.  The  results  are  summarized  in  Table  III-l  which  gives  the  pulse 
delay  time,  width  intensity,  and  area  under  the  pulse. 

Figure  in-4  (a,  b,  c)  shows  the  parametric  fit  of  the  delay  time  xq  as  a  function  of 
the  attenuation  coefficient  p.  to  the  quadratic  function 

tjj=(a+Pp  +  Yli^)xo  .  (19) 


The  curve  fit  is  for  the  calculated  points  that  result  in  the  curves  in  Fig.  III-2. 

Table  III-2  gives  the  corresponding  parameters  that  define  the  quadratic  fits. 

It  has  been  claimed  (Refs.  1,  11)  that  a  finite  attenuation  coefficient  affects  the  ^ 

development  of  SF  pulses  by  restricting  the  active  cooperating  volume  by  limiting  the 

cooperation  length  to  •^.  We  tested  this  assumption  by  comparing  the  delay  time  xdmc 

obtained  from  the  Monte  Carlo  solution  of  the  Maxwell-Bloch  equations  with  the  result,  # 

N 

designated  xdv  for  this  purpose,  obtained  from  Eq.  (17)  assuming  N  ->  ,  or  in 


t  X 

normalized  units  of -= — andy  , 
^SF  1 


‘’t>v~ 


=  a 


1  +  bln 


(sMf)-  - 


(20)^ 


At  p.  =  1,  Xdv  from  equation  (17)  for  different  N  was  normalized  to  the  Monte  Carlo 
result,  Xdmc  (shown  in  Fig.  111-2). 


1 


XDV  is  in  units  of  tsF<  and  the  factor  p  on  the  extreme  right  comes  from  the  decrease  of  the 
cooperation  number  used  in  the  definition  of  tsp. 
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Figure  111-2.  SF  Delay  time  to  as  a  function  of  linear  attenuation  coefficient  (i. 
Points  obtained  from  Monte  Carlo  calculation. 
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Figure  Ill'S.  Normalized  area  under  SF  pulsea  obtained  from 
Monte  Carlo  calculation  as  a  function  of  p. 
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Table  III-1.  SF  Pulse  Ctiaracteristic  Parameters 


Attenuation 
N  Coefficient  m 


20 

77 

40 

155 

80 

338 

Pulse 

Intensity  Width 


17.5  0.0329 


20.0  0.0194 


.00439 


67.0  0.00108 


125.0  0.000198  70.0 


240.0  0.0000385 


23  0.0247 


0.0168 


0.00381 


0.00104 


.000211  75. 


0.0000231  180.0 


0.0189 


Area 

Under  Pulse 


.378 


0.272 


.079 


.038 


0.014 


0 

0.0000385 

156.0 

■ 

.006 

■ 

0.333 


0.235 


0.08 


0.036 


.016 


.00416 


0.302 


Normalized 

Area 


1.0 


0.72 


0.21 


0.10 


0.04 


.016 


1. 


0.24 


0.11 


0.048 


0.012 


1.0 


2 

37 

0.0126 

16.5 

0.208 

0.69 

10 

58 

0.00353 

42.0 

0.085 

0.28 

20 

110 

0.000991 

33.0 

0.033 

0.11 

40 

190 

0.000181 

64.0 

0.0116 

0.038 

0.000016  i  232.0 


0.0037 


0.012 


II 


UNEAR  ATTENUATION  COEFFICIENT,  \i 

Figure  lll-4a.  Parametric  fit  to  data  ehown  In  Fig.  IV*2  for  N  s  10^. 
The  best  fit  Is  given  by  Xo=  (1278+220|i+1/4Ki»)  1.17  x  10"*. 


Figure  lll-db.  Parametric  fit  to  data  shown  In  Fig.  IV-2  for  N  s  10^. 
The  best  fit  is  given  by  Xo®  (307+37ji+i/4n2)  6.93  x  itr*. 
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UNEAR  ATTENUATION  COEFFICIENT,  \L 

11-23-90-10  ^ 

PIgura  lll^.  Parametric  fit  to  data  ahown  In  Fig.  IV*2  for  N  s  io^. 
The  best  fit  la  given  by  xd  ^  9.1  xl0~2  (327<i>35|i  1/4 


Table  Ilh2.  Parameters  for  Delay  Time  as  a  Function  of  Attenuation 


N 

P 

7 

to 

102 

1278.00 

220.20 

0.25 

1.17x10-2 

IqS 

306.89 

37.25 

0.25 

6.93x10-2 

10^ 

326.85 

34.87 

0.25 

9.11  xlO-2 

► 

The  comparisons  are  shown  in  Hg.  ni-5(a-c)  with  the  heavy  line  representing 
the  Monte  Carlo  calculation  Tdmc  and  the  dashed  line  the  reduced  volume  result  Xdv- 
^  For  N  =  10^  [Fig.  ]ll-5(a)],  the  two  curves  intersect  at  p  ~  45  and  Xov  approaches  zero 

since  approaches  1.  Gearly  the  meaningful  range  of  p.  is  p.  <  20  for  N  =  10^.  For 

hi^er  N  values  this  restriction  occurs  at  hi^er  p  values.  The  more  interesting  result  of 
this  (xxnparison  is  that  the  volume  restriction  by  p  which  is  the  basis  for  the  calculation  of 
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Figura  lll-5(a-f).  A  plot  of  SF  delay  tlmea  tpv  ^OMC  ••  ■  function  of  the  linear  attenuation  coefficient  p.  The  delay  time  tpy 
(dashed  line)  was  obtained  from  the  reduced  volume  assumption  [Eq.  (2)1  and  tomc 

simulations.  The  cooperation  number  N  Is  10*  In  (a),  10®  In  (b),  10*  In  (c),  10*  In  (d),  10*  In  (e),  10^  In  (F). 


gives  a  larger  value  for  the  delay  time.  If  we  assume  that  the  Monte  Carlo  calculation  is 
a  better  representation  of  the  physical  process  and  therefore  represents  the  conect  value,  the 
volume  reduction  or  restriction  result  Xdv  has  to  be  modified.  From  the  three  figures 
showing  the  comparison  between  Xdv  tuid  Xdmc>  see  that  in  the  p.  <  20  range  by 
reducing  the  attenuation  coefficient  p  by  a  factor  of  3  an  approximate  value  for  xdv  closer  to 
Xi^c  can  be  obtained  from  Eq.  18.  The  factor  increases  for  higher  p  values,  for  example: 
at  N  s  1(P  and  p  ~80  the  required  factor  to  make  the  correction  would  be  closer  to  4. 
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IV.  SUPERFLUORESCENCE  IN  <50Co 


f 

^  A.  CHARACTERISTIC  PARAMETERS 

In  diis  section  we  apply  the  results  of  our  earlier  studies  to  examine  the  possibility 
of  SF  in  ®^Co.  We  first  determine  the  SF  time  xsf  for  ^Co.  We  use  this  to  generate 
^  delay  time  xd  values  (obtained  from  different  theories)  and  we  then  check  the  SF  condition 

^XspXp  <  to  determine  the  bound  on  the  only  free  (unknown)  parameter  a, 

characterizing  the  inhomogeneous  broadening.^  Finally  we  check  our  preliminary  results 
with  the  Maxwell-Bloch  Mcmte  Carlo  calculations.  The  characteristic  parameters  are  givoi 
r  in  Table  IV-1,  and  the  partial  energy  diagram  for  is  shown  in  Hg.  IV-l 


¥ 


1  is  the  dqthasing  liinc  due  to  inhomogeneous  broadening,  etc. 

IV-1 

■ _ 


Table  IV*i.  Characteristic  Parameters  for 


Parameter 

Value 

Symbol 

Transition  energy 

58.6  ke  V 

EO 

Wave  length  on  resonance 

02A{2x10-®cm) 

X 

Natural  Wetime 

628  sec. 

To 

ReooiHess  fraction 

0.304 

f 

Lirtear  attenuation  coefficient 

12  cm"^ 

Intemai  conversion  coefficient 

48.3 

a 

Particte  density  (in  solid  form) 

8.97  X  10^2  cm”3 

P 

Thermal  neutron  cross-section  ^  60^^ 

20  bams 

On 

4+ 


277.20 


10.47  m 

0.24  %  5.2714  y 

100% 


5-1 0-91 -1M 

Figure  IV-1.  Level  Structure  of  ^*^0027 

The  characteristic  SF  time,  assuming  an  effective  length  limited  by  the  attenuation 
coefficient  and  ^proximately  4^1  (see  Chapter  IE),  is  given  by 

Sjtt  8jt(l+a)x 

— 2 - ’ 

3fpX  /  3fpx  (4/p) 

which  for  gives  tsf  =  7.1  sec.  In  units  of  xsF  lifetime  is  88  and  the  natural  decay 
rate  F  =  1.1  x  10-2 

The  delay  time  has  been  derived  in  different  developments  with  different 
dependences  on  the  number  of  cooperating  nuclei  N.  The  number  of  cooperating  nuclei  is 
given  by 

N  =  pX2/  =  9  X  104, 

using  a  value  of  0.25  cm  for  /.  For  the  delay  time,  Bonifacio  and  Lugiato  and  Gross  and 
Harocheuse 

Xj,  =  ixsF/n(N)  ;  (21) 

on  the  other  hand.  Polder,  Schuurmans,  and  Viehen  derived  the  expression 


1 

V  =  A»(27cN)^ 

D  4  SF  »- 
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for  that  quantity.  They  claim  that  the  /n  ^(N)  dependence  is  due  to  taking  into  account 
prc^agadon  effects  that  Bonifacio  and  Lugiato  do  not  include  in  their  model.  From  our 
Monte  Carlo  calculations  in  studies  of  the  MaxweU-Bloch  equations,  including  random 
sources  and  {mogadon  effects,  we  found  that 

Tjj=a[l+b/n(N)]/n(N)  (23) 

with  a  a:  2.43  and  b  =  0.02.  These  three  results  are  compared  in  Fig.  rV-2. 


NUMBER  OF  RADIATORS,  N 

5a-91>2Ma 

Figure  IV-2.  Comparison  of  Pulse  Delay  Times  calculated  from  different 
theories  as  in  (a)  Bonifacio  and  Lugiato  (Ref.  9)  and  Gross  and  Haroche 
(Ref.  10)  results  (tp  of  equation  8),  (b)  Polder,  Schuurmans  and  Vrehen 
(Ref.  12)  result  ()^  of  equation  9),  (c)  Haake  and  Relbold  (Ref.  13) 

result  (XqOI  equation  10). 


For  Large  N  the  effect  of  the  In  2(N)  term  becomes  important  but  for  N  =  105  this 
term  contributes  only  about  20  percent  of  the  In  (N)  term's  contribudon.  Table  rV-2  gives 

die  values  of  := — ,  and  ^  for  different  values  of  N. 


-SF  "SF 


‘SF 
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Table  IV>2.  Comparison  of  Pulse  Delay  Times  calculated 
as  a  function  of  N 


N 

/n(N) 

V 

V 

it 

^SF 

102 

4.6 

2.3 

12.2 

2.6 

103 

6.9 

3.4 

19.1 

4.8 

10^ 

9.2 

4.6 

26.5 

mam 

105 

11.5 

5.8 

34.4 

11.1 

10® 

13.8 

6.9 

42.8 

15.3 

107 

16.1 

8.0 

51.8 

20.1 

10® 

18.4 

9.2 

61.2 

25.6 

B.  BOUNDING  CALCULATION 

Using  tsF  =  7.1  sec,  wc  generate  Xd  and  tp  for  ^Co  for  different  values  of  N 
shown  in  Table  IV-3. 


Table  IV*3.  Delay  Times  Calculated  two  different  ways 


N 

tD(sec) 

^(sec) 

102 

16.3 

86.6 

10® 

24.1 

135.6 

10^ 

32.7 

188.2 

10® 

40.8 

244.2 

10® 

49.0 

303.9 

107 

56.8 

367.8 

10® 

65.3 

434.5 

All  of  the  calculated  delay  times  are  smaller  than  the  natural  lifetime  (628  sec),  and 
so  in  the  absence  of  attenuation  and  dephasing  (inhomogeneous  broadening)  an  inverted 
population  of  should  produce  SF  pulses.  We  still  have  to  include  these  effects  for  a 
nxne  realistic  estimate  and  this  will  be  done  in  the  next  section. 
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A  relationship  exists  that  allows  us  to  put  an  upper  bound  on  the  inhomogeneous 
broadening  parameter  in  terms  of  p  and  a,  which  are  known  (Table  FV-l).  In  the 
approximation  using  the  Xj)  derived  by  means  of  semiclassical  or  mean  field  theory  for  the 
delay  time  we  have  [fnnn  Eq.  (16),  Chapter  II] 


(j+l)(a+l)(a+l)< 


4jtV2/n(N)  ’ 


which  with  a  =  48.3  and  p  =  12  and  N  =  lO^  gives  a  <  27  (for  N  =  lO^,  a  <  35)  as  the 
requirement  on  inhomogeneous  broadening  in  order  to  achieve  SF.  When  we  use  in  the 
estimation  the  requirement  becomes  a  <  10  for  N  =  10^  and  a  <  14  for  10^.  A  similar 
calculation  for  ASE  in  ^Co  gives  a  <  1.4  as  the  requirement  on  inhomogeneous 
broadening. 


In  the  directions  of  the  crystal  where  the  Brarmann  effect  reduces  the  attenuation  to 
p  »  0  the  requirement  for  SF  is  a  <  1 12  for  N  =  10^  when  xd  is  used  in  the  estimate.  All  of 
these  estimates  depend  on  full  inversion  at  t  =  0.  To  include  finite  inversion  times  and  treat 
all  of  the  pertinent  parameters  in  a  consistent  way  we  require  the  solution  of  the  Maxwell- 
Bloch  equations  as  described  in  Refs.  4  and  8.  This  will  be  discussed  in  the  next  section. 
As  a  final  summaiy  of  the  bounding  study,  we  present  curves  that  bound  the  values  p  and  a 
for  different  cooperation  numbers  N. 


Figures  IV-3(a,  b)  give  the  SF  conditions  according  to  bounds  determined  by 
Eq.  (16).  For  each  curve  representing  a  particular  value  of  nuclear  density  p  the  area  below 
and  to  the  right  gives  the  acceptable  values  of  a  and  p.  Figure  IV-4  gives  the 
corresponding  results  for  ASE  according  to  bounds  determined  by  Eq.  (IS).  Figure  rV-4 
for  ASE  is  drawn  to  the  same  scale  as  the  SF  curves  in  IV-3  (a,  b).  Note  the  restricted 
parameter  space  for  ASE  as  compared  to  SF. 


C.  MAXWELL-BLOCH  RESULTS 

TTie  interplay  of  the  important  physical  parameters  in  dynamical  full-time  dependent 
treatment  of  the  physical  phenomena  is  needed  fca*  a  realistic  assessment  of  the  feasibility  of 
SF.  In  this  section  we  present  our  results  on  the  SF  pulse  emission  in  ^^Co  from  the 
Maxwell-Bloch  Monte  Carlo  calculations.  Using  parameters  given  in  Table  IV-l  as  the 
nominal  values,  we  examined  the  possibility  of  SF  by  studying  the  pulse  as  a  function  of 
the  inhomogeneous  broadening  parameter  a,  and  attenuation  coefficient  p.  Some  SF  pulses 
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INHOMOGENOUSniMOENINGPMMlinBI  a 

II  23'«>'22 

Figure  IV-4.  A  plot  of  ^  versus  a  for  ASE  at  the  boundary.  Each  curve  Is 
calculated  tor  a  particular  value  of  the  density  p  as  martced.  All  the 
points  in  the  region  to  the  left  of  the  curve  give  the  acceptable 
values  of  a  and  |x  for  ASE. 


calculated  for  ^Co  with  the  Maxwell-Bloch  equations  described  in  Ref.  8  are  shown 
in  Fig.  IV-5,  IV-6  and  rV-7.  Pulse  peak  intensity  I©,  delay  time  to  (in  seconds),  and 
the  integrated  pulse  intensity  per  nucleus,  for  the  different  values  of  parameters 
the  attenuation  coefficient  (in  cm~l),  and  a,  the  inhomogeneous  broadening  parameter, 
are  given  in  Table  IV-4.  For  these  calculations  we  assumed  1  =  1  cm,  giving 
N  =  pX^l  =  1.1  X  105  and,  from  Eq.  (2),  xsf  =  2.4  seconds.  Optimizing  for  the 
geometrical  divergence  and  the  diffraction  limits  on  a  cylindrical  active  region  of 
length  1  and  diameter  d,  requires  d=//X.  For  1  =  1  cm  and  A,  =  2.0  x  10“9cm,  we  get 
ds  4.5  X  10-5  cin  and  a  beam  divergence  of  6d  =  9  x  10-5  radians. 

The  possibility  of  observing  nuclear  SF  is  considered  by  comparing  the  SF 
emission  as  characterized  in  Table  rV-4  with  the  spontaneous  emission  from  the  natural 
decay  of  ^Co.  The  two  experimental  techniques  considered,  A  and  B,  are  shown  in 
Fig.  IV-8.  In  experiment  A,  a  sample  composed  of  small  single  crystals  of  ^Co, 
randomly  oriented,  is  placed  in  a  thermal  neutron  reactor  and  the  population  inverted.  The 
sample  is  then  placed  in  a  47C  detector  and  the  count  as  a  function  of  time  measured. 
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I.  SF  pulses  calculated  with  the  Maxwell-Bloch  equations  In  a  Monte  Carlo  simulation  of  the 

vacuum  fluctuations.  The  parameters  arey  *  10,  r2  *  10“®,  Fs  ss  3.75  x  10"* 
and  N,  \i  and  a  as  Indicated  In  the  figure. 


SF  pulses  calculated  with  the  Maxwell-Bloch  equations  In  a  Monte  Carlo  simulation  of  the 

vacuum  fluctuations.  The  parameters  are  y  s  10,  r2  s  10~*,  Fs  s  3.75  x  10”® 
and  N,  p  and  a  as  Indicated  In  the  figure. 


Table  IV-4.  SF  Pulse  Parameters  calculated  for  soco  from  the  Maxwell>Bloch 
equations.  Shown  are  the  SF  pulse  delay  time  tq.  normalized  pulse 
Intensity  (photons/second)/N  In  parenthesis  (  )  and  imegrated 
normalized  Intensity  (Photons/^)  |n  brackets  [  ]  for  different 
values  of  the  attenuation  p  (cm-f)  and  inhomogeneous 
broadening  a.  assuming  N  s  1.1  x  lOS  cooperating  nuclei 
and  a  pumping  rate  of  4.2  s-'*. 


N.  a 

P 

0 

5 

10 

20 

50 

0 

119 

(0.018) 

[0.348] 

136 

(0.076 

[0.163] 

167 

(2x10-^) 

[0.0703] 

122 

(5.47x10-5) 

[2.17x10-3] 

48 

(5.96  X  10-2) 
[1.0x10-5] 

6 

190 

(0.0055) 

[0.118] 

221 

(1  X  lO-^) 
[0.0376] 

252 

(4.1  X  10-5) 
[2.75  X  10-3) 

133 

(4.77x10-7) 

[1.95x10-5] 

48 

(5.16x10-5) 

[8.62x10-*] 

12 

269 

(0.0017) 

[0.0417] 

381 

(3.85x10-5) 
[3.38  X  10-3] 

252 

(4.41  xlO-7) 
[3.18  X  10-5] 

129 

(5.08x10-5) 

[2.07x10-7] 

3 

(5.89x10^10) 
[1.02  X  10-®] 

Figure  IV-8.  Two  experimental  setups  for  detecting  nuclear  SF.  In  experiment 
A,  a  polycrystalline  sample  is  used  and  both  SF  and  NSO  give  isotropic 
radiation  distributions;  while  In  experiment  B  an  acicular  geometry 
provides  a  highly  directional  SF  beam. 

EXPERIMENT  A  EXPERIMENT  B 
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No  atteiiq>t  is  made  at  angular  discrimination  of  the  SF  emissitm  from  the  natural  decay  and 
only  the  measurement  of  the  emission  rate  is  used  to  detect  SF.  In  experiment  B  the  sample 
is  prepared  in  a  cylindrical  geometry,  irradiated  and  placed  in  front  of  a  collimator  so  that 
the  directional  beam  characteristics  can  be  used  to  discriminate  against  the  natural  decay. 
The  ratio  of  the  photon  count  from  SF  pulses  divided  by  the  spontaneously  emitted  photons 
in  the  natural  decay  process  integrated  over  all  angles  (detected  in  a  4tc  detector)  is  2.9 
when  a  s  0,  p.  =  0  (case  1);  1.0  x  10-^  when  a  =  50  and  p,  =  0  (case  2);  and  8.6  x  10-6 
when  a  =  20,  p  =  12  (case  3).  We  also  lowered  the  pumping  rate  to  y  =  10~6  (equivalent  to 
2.1  X  10^1'^  neutrons/cm^)  and  repeated  the  a  =  0,  p.  =  0  calculation  (case  4)  with  the 
result  1.6  X  10-8.  if  advantage  is  taken  of  the  directionality  of  the  SF  pulse  in  the 


experimental  design  then  the  small  divergence  of  the  SF  beam  compared  to  the  isotropic 
distribution  of  the  natural  decay  would  permit  discrimination  and  reduction  of  the  effect  of 
the  natural  decay,  in  the  axial  direction  as  compared  to  the  SF  pulse,  by  a  factor  of 


^  =  2.0  X  10  The  photon  count  from  the  SF  pulse  compared  to  the  natural  decay 


would  increase  to  1.5  x  10®  for  a  =  0,  p  =  0  (case  1),  4.9  x  lO^  for  a  =  50,  p.  =  0  (case  2), 
and  4.2  x  10^  for  a  =  20  and  |X  «  12  (case  3).  For  case  4  we  obtained  25  as  the  ratio  of  SF 


to  natural  photons.  These  results  are  tabulated  in  Table  IV>5.  In  suinmary,  for  a  Borrmann  • 

direction  along  the  propagation  axis  (case  1),  SF  could  be  detected  on  the  basis  of  time 

development  of  the  decay,  and  no  discrimination  based  on  angular  distribution  is 

necessary.  In  case  2,  again  along  the  Borrmaiui  direction  but  now  with  inhomogeneous 

broadening  of  a  =  50,  angular  discrimination  would  have  to  be  used  to  detect  the  SF.  In  • 

case  3,  without  the  Borrmann  effect  present  and  with  inhomogeneous  broadening  of 

a  =  20,  the  SF  could  also  be  observed  with  the  help  of  angular  discrimination  obtained 

through  collimation.  In  case  4,  with  the  reduced  pumping  rate,  both  the  Borrmann  effect 

and  reduced  inhomogenous  broadening  are  required  for  the  observation  of  SF.  • 


Figures  IV-9  and  IV- 10  show  calculated  pulse  shapes  in  the  presence  of  natiual 
decay  for  case  1  without  assuming  angular  discrimination,  and  case  2  with  angular 
discrimination  of  only  6  =  10-^  radians  more  easily  available  experimentally  instead  of  the 
lower  limit  permitted  by  the  divergence  of  the  SF  beam  of  ©d  =  9  x  10-5.  jn  both  figures, 
(a)  shows  separately  the  photon  intensity  calculated  for  both  processes  and  (b)  the  sum  of 
the  intensity  of  the  two  processes. 


IV.12 


Table  IV*5.  Predictions  for  *^Co  Experiments 

Total  Count  Ratio 
Nsf/Nnso 


Experiment  A 

Experiment  B 

poiycrystal 

whisker 

single  crystal 

Photons 


6-4-91 -4m 


t  (seconds) 


t  (seconds) 

6-4-91 -4am 


Figure  IV-9.  SF  and  Natural  Spontaneoua  Decay  (NSD)  In  Assuming 
no  attenuation  (u  s  0),  no  Inhomogeneous  broadening  (a  s  0),  and 
no  Directional  Discrimination.  The  scale  Is  linear  for  both  axes. 

The  Intensity  Is  given  In  photons  per  second  per  emitter. 
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We  now  have  to  project  the  results  of  our  calculations  onto  the  real  world 
situation  and  consider  the  experimental  realization  of  our  predictions  for  observing  SF  in 
^^o.  The  calculations  leading  to  the  results  shown  in  Table  rV-4  assumed  a  pumping 
rate  in  normalized  units  of  Y=  ^0.  This  corresponds  to  a  thermal  neutron  flux  of 
J  =  2.1  X  1023  cm~2s-l  which  can  certainly  be  obtained  as  a  result  of  nuclear  explosions 
(Ref.  14),  but  is  out  of  reach  of  available  continuously  operating  reactors.  Table  rV-6 
shows  the  calculated  degradation  of  the  SF  pulse  intensity  as  the  pumping  rate  decreases. 
The  lower  pumping  rates  with  y  =  10^.  10“^  require  fluxes  that  were  considered  within 
reach  of  practical  reactor  designs  in  1975,  and  so  should  be  at  least  physically  attainable 
from  the  point  of  view  of  reactor  designers. 

Table  IV>6.  Total  Integrated  Intensity  of  SF  pulses  as 
a  function  of  pumping  rate. 


The  Borrmann  effect  has  been  observed  and  predictions  indicate  that  suppressitm  of 
attenuation  by  two  orders  of  magnitude  is  certainly  possible  under  the  right  conditions 
(Ref.  15).  Inhomogeneous  broadening  has  been  discussed  in  the  literature,  and  theoretical 
predictions  using  a  static  nuclear  dipole-dipole  interaction  model  give  a  value  of  a  =  10^  for 
this  isomer  (Ref.  16).  Subsequent  measurements  by  three  different  groups  (Refs.  17-19) 
of  self  attenuation  in  a  lO^Ag  source  gave  a  »  30-50.  In  addition  to  this,  relaxation  was 
proposed  as  a  possible  mechanism  (Ref.  19)  to  explain  the  discrepency  between  these 
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Pumping  Rate  I  Neutron  Flux 


Y(Vtsf)  J(crn-2-sri) 


0.002 


2.1  X  1023 


2.1  X  1022 


2.1  X  102^ 


2.1  X  1020 


4.2x10^2 


2.1  X  lO’S 


2.1  X  10^® 


2.1  X  10^7 


Integrated  lntensity/|^ 


Case  1 
u-0,a-0 

Case  2 
p>0,a«50 

Case  3 

U«  12,a-20 

0.348 

1.1  X  10-8 

2.0  X  10-7 

0.348 

6.5  X  10-6 

6.7  X  10-8 

0.264 

8.3x10-8 

5.0x10-8 

0.146 

3.1  X  10-9 

7.2x10-9 

8.95  X  10-8 

7.1  X  10-^0 

1.4x10-9 

5.1  X  10“^ 
13  X 10-7 
1.6x10-8 


measurements  and  the  calculated  value  of  10^  from  the  static  dipole-dipole  interaction 
(Ref.  16).  We  discuss  the  effect  of  relaxation  on  SF  in  another  publication  (Ref.  8). 

D.  CONCLUSIONS 

This  paper  considers  a  generalized  version  of  the  Haake-Reibold  model  of 
superfluorescence  in  a  multilevel  system  of  emitters,  having  up  to  four  energy  levels.  The 
generalized  version  takes  into  account  finite  pumping  rates,  spatial  attenuation  of  the 
electromagnetic  field,  dephasing,  and  population  depledcm  due  to  processes  such  as  internal 
conversion  and  spontaneous,  but  noncollective,  transverse  emission. 

In  this  paper  we  study  the  dependence  of  emitted  pulse  shapes  on  the  statistics  and 
time  dependence  of  the  noise  sources.  Monte  Cailo  results  obtained  by  averaging  25  to 
100  pulses,  which  seem  to  be  adequate  samples,  also  depend  upon  the  noise  source 
vaiiance. 

We  find  diat  the  spatial  attenuation  parameter  |i  reduces  the  pulse  peak  and  shifts  it 
to  longer  times.  The  peak  will  decrease  gradually  but  will  not  disappear  if  dephasing 
mechanisms  that  reduce  the  macroscopic  dipole,  built  up  by  the  noise  source,  do  not  exist 
and  population  depletion  mechanisms  do  not  destroy  the  inversion. 

Small  values  of  the  dephasing  factor  have  no  effect  on  the  emitted  pulse  delay  time 
or  intensity,  but  when  the  dephasing  rate  due  to  either  homogeneous  or  inhomogeneous 
broadening  is  con^aiable  to  the  superfluorescence  decay  rate,  l/xsF>  the  delay  time  starts  to 
increase  with  increasing  values  of  the  dephasing  factor  until  it  reaches  a  maximum,  after 
which  it  decreases  rapidly.  Consistently,  the  inhomogeneous  effect  is  stronger  than  the 
homogeneous  effect  for  tiie  same  amount  of  broadening. 

Because  of  spatial  attenuation,  for  maximum  intensity  of  the  pulse  radiated  out  of 
the  active  region,  the  choice  of  the  region's  length  should  be  governed  by  the  values  of  the 
attenuation  factor  ^  and  the  coupling  constant  g]  associated  with  coupling  between  the 
emitters  and  the  electromagnetic  field.  The  optimum  length  will  depend  on  details  of  the 
propagating  electromagnetic  field  inside  the  active  region. 

The  theory  developed  in  this  paper  is  applied  to  the  58  keV  transition  in  ®0co.  It  is 
found  that  for  sufficiently  intense  pumping  with  thermal  neutrons  an  inversion  in  can 
be  achieved  and  a  SF  pulse  detected  experimentally. 
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APPENDIX  A 

THE  TRAMMELL-HANNON  MODEL  OF  AMPLIFIED 
SPONTANEOUS  EMISSION  (ASE) 


In  this  appendix  we  discuss  pulse  propagation  through  an  amplification  medium, 
calculate  pulse  shapes  as  a  function  of  gain  and  attenuation,  and  consider  the  appropriate 
form  of  the  threshold  conditions  for  pulse  amplification.  Our  calculations  are  based  on  the 
Trammell  and  Hannon  (Ref.  A-1)  model  of  ASE  which  has  certain  severe  limitations.  In 
this  appendix  we  modify  the  model  to  remove  the  small  signal  restriction.  The  study  of 
ASE  is  important  because  it  bears  on  the  question  of  whether  ASE  or  SF  pulses  should  be 
expected  from  an  inverted  medium  and  therefore  on  the  selection  of  candidate  nuclei  for  a 
y-ray  laser. 

Trammell  and  Hannon  (Ref.  A-1)  doived  an  expression  for  a  pulse  propagating 
through  a  gain  medium  of  length  /  and  attenuation  coefficient  |X.  Assuming  an  incoming 

pulse  of  unit  amplitude  and  line  width  F  and  corresponding  lifetime  ^  they  concluded 
that  the  outgoing  pulse,  after  passing  through  the  medium,  would  have  the  form 

Kt)  =  I  A  (/.  t)  I  ^  =  e-'^  I  rt  An* 

where  Ig  is  the  modified  Bessel  function  of  order  zero  and 

Kq  =  2n}^  fg  (1+a)'^  (1+a)'^  .  (A-2) 


X 

With  X=-^,  and  X  the  wavelength  corresponding  to  the  energy  on  resonance; 

f  is  the  recoilless  fraction  of  the  transition; 
g  is  the  coupling  of  the  electromagnetic  field  to  the  nucleus; 


a  is  the  internal  conversion  coefficient; 
a  is  the  inhonrageneous  broadening  parameter. 


An*(t) 


is  the  inversion  at  time  t; 


ni(t),  n2(t)  are  the  populations  and  g2,  gi  are  the  statistical  factors  for  the 
two  states. 
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Trammell  and  Hannon  assumed  that  for  nuclear  transitions  of  interest  to  the  y-ray 
laser  problem  it  was  valid  to  replace  An*(t)  by  the  function  njCO)  [2e  —  1].  This 

assun:q)tion  means  that  the  change  in  the  level  populations  is  determined  principally  by  the 
natural  lifetime  and  is  strictly  true  only  in  the  limit  of  small  gain  Kq  1.  If  the  gain  is  large 

the  approximation  breaks  down. 

It  is  possible  to  deal  with  cases  in  which  the  gain  is  larger  by  replacing  (1)  with 
solutions  of  the  nonlinear  Haake-Reibold  equations,  which  involve  both  the  electric  field 
and  the  polarization  and  which  were  solved  numerically  in  Ref.  A-2.  However,  a  simpler, 
phenomenologically  motivated  numerical  procedure,  which  operates  directly  on  the 
intensity  in  the  form  given  by  (1),  can  be  used  to  estimate  the  behavior  of  pulses  as  the  gain 
is  increased. 


The  algorithm  for  this  purpose  is  based  on  the  following  considerations.  The 
iiutially  inverted  population  must  be  equal  to  the  total  number  of  photons  emitted  over  time, 
i.e.. 


N. 


=/ 

o 


A  n  dt 


(A-3) 


Initially 


ni(0)  =  No.n2(0)  =  0. 


I 

(A-4) 


At  any  time  t,  after  a  time  increment  At  the  corresponding  decrease  An  (t)  in  the  populatitm 
is  the  result  of  two  processes:  (1)  the  noninter^tive  emission 


nCOoTAt  , 

which  is  due  to  the  natural  decay  of  individual  emitters;  (2)  the  interactive  emission,  which, 
accOTding  to  (A-1)  if  absorption  is  ignored  (p.  =  0)  is  given  by 


where 


X(/,t-At)At  , 

X(/,t)  =  e-^I^[^^/  rtAn*(t)]  . 


I 


(A-5) 


The  algorithm  begins  by  setting  the  initial  conditions  Eq.  (A-4).  Then  it  goes 
through  the  following  steps  to  upgrade  the  quantity  X  (/,  t)  whenever  the  time  t  is 
incremoited  by  an  amount  At: 
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1 .  An  =  no(t)  FAt  +  X  (/,  t  -  At)  At; 

2.  ni  (t)  =  ni  (t  -  At)  -  An  (t); 

3.  n2  (t)  =  n2  (t  -  At)  +  An  (t); 

4.  An*  (t)  =  ni  (t)  -  n2  (t). 

5 .  Substitute  the  result  of  step  4  into  Eq.  (AS)  above. 

Then  die  ASE,  taking  into  account  absorption,  is  given  by 

“  |A(/.t)|2  =  X(/.t)e-W. 

which  follows  from  Eqs.  (A-1)  and  (A-S). 

For  the  nuclear  parameters  shown  in  Table  A-1  gain  is  small  and  the  pulse  shown  in 

Fig.  A-1  for  g  =  1  appears  close  to  the  value  *  *  as  determined  by  Trammell  and 

Hannon  but  the  peak  is  not  exactly  at  the  value.  We  have  found  that  for  higher  g  (or  larger 
wavelength  X)  the  pulses  become  asymmetric  and  shift  toward  shorter  times,  as  shown  in 
Fig.  A-1.  The  other  parameters  used  to  calculate  the  pulses  shown  in  Fig.  A-1  are  given  in 
Table  A-L 


Table  A*l.  Nominal  Parameters  used  In  the  Calculations 


For  a  constant  g  and  increasing  nuclear  density  p  the  pulses  are  shown  in  Fig.  A-2. 
Note  that  the  pulses  disappear  between  p  =  10^  and  10^^.  This  corresponds  to  the  region 
KoAn(o)  =  13.15  to  1.315.  The  pulse  should  disappear  even  for  p  =  0  when  the  condition 

I^(t)(^Ko/  rtAn*(t))er‘<  1  is  met.  For  this  condition  any  effect  of  stimulation  is 
overcome  by  the  natural  decay. 

We  also  investigated  the  decrease  in  ASE  as  a  function  of  p.  Figure  A-3  shows 
pulses  for  p  =  ICp  and  p  =  10^  and  different  values  of  p.  Trammell  and  Hannon  derived 
a  modiHed  Schawlow-Townes  condition  for  the  existence  of  ASE  pulses.  Their 
nndification  is  to  aiiqrlify  the  effect  of  p  by  writing  the  Schawlow-Townes  ccxidition  as 

K  =  KoAn*-Mp>0  .  (A-6) 


The  factor  M  in  front  of  p  results  from  the  decrease  in  inversion  in  pulse  lasers 
during  the  pulse  emission.  When  subject  to  the  assumption  made  by  Trammell  and 
Harmon  that  An*(t)  =  ni(0)(2e^  -  1),  it  was  determined  that  M  =  21.  We  do  rrot  make  this 
assumption  in  our  simulation  studies.  In  our  case  An*(t)  is  calculated  including  the 
reduction  by  the  ASE  tnechaiusm  as  well  as  natural  decay.  We  can  determine  the  threshold 
condition  for  lasing  by  determining  the  effect  of  the  factor  e~^  cm  the  amplification  directly 


Ending  the  value  of  p  for  which 


c*^lj(yKo/  rt  An*(t))]  e-*^ 


becomes  less  than  one. 


We  first  calculate  I(t)  from  equation  (A-1)  for  p  =  0  and  then  determine  the  value  of  p 
needed  to  reduce  I(t)  at  maximum  to  1.  We  compare  this  to  the  p  calculated  from  the 
trxxlified  Schawlow-Townes  ccmdition  with  M  =  21  [Eq.  (A-6)]. 


Table  A-2  gives  the  results  with  colurrm  1  giving  the  peak  intensity  I(t)  for  p  =  0, 
colurrm  2  giving  the  threshold  value  of  p  calculated  from  Eq.  (A-6),  column  3  giving  the 
value  of  p  required  in  Eq.  (A- 1)  to  reduce  the  value  in  colurrm  1  (in  Table  A-2)  to  1,  and 
colurrm  4  gives  the  value  of  p  required  to  produce  threshold  condition  in  Eq.  (A-6)  with  p 
obtained  from  colurrm  3.  We  use  the  symbol  p  for  the  attenuation  coefficient  obtained  from 
the  modified  Schawlow-Townes  condition  Eq.  (A-6)  and  p  for  the  value  required  in  the 
exponent  in  Eq.  (A-1)  to  reduce  the  pulse  intensity  to  1. 
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Figure  A-2a.  Logarithmic  plot  of  the  ASE  pulses  for  g  s  1  and  different  values 
of  the  Initial  Inversion  density,  p  s  10*®,  10*’’ ,  10**,  10**. 


TIME 


Figure  A>2b.  Plot  of  the  ASE  pulses  for  values  of  the  Initial  Inversion  density 
In  the  threshold  region  between  psio^®  and  p  sio*®  cm~*. 


Table  A>2.  Threshold  conditions  on  ASE 


_ m _ f21 _ . _ Idl 


Critical  attenuation  coefficient  p  (ctn*^) 
calculated  using  two  approaches 

Factor  multiplying  p  in 
equation  3  for  no  gain 

density 

p(cnr3; 

KoP 

peak  intensity** 
(sec~^)whenu-0 

From  Modfied  S-T 
Condition**  (p) 

M 

<0 

CM 

O 

1.315  X  10^ 

5x10^ 

650 

54.5 

250 

1022 

1.315  X  103 

10® 

65 

20.7 

68 

1.315  X  102 

10^ 

7.5 

4.6 

32 

IE3 

1.315  X  10 

1.7 

0.8 

0.53 

21 

PmK  intensity  obtained  from  Figure  A-2. 
Critical  m  calculated  from  Eq.  (A-6). 


***  Critical  m  calculated  from  Eq.  (A-1 )  assuming  no  amplification  |A(/  , !)(  =  1  • 

From  Table  A-2,  we  note  that  if  we  want  to  use  the  Schawlow-Townes  condition 
the  factor  multiplying  p  in  Eq.  (A-3)  has  to  increase  for  systems  with  high  gain  or  high 
density  to  about  250  for  p  =  10^^  cm~^.  Thus  for  ASE,  M  is  a  function  of  K^An*,  but  not 

a  linear  function.  For  low  gain  M  >=  21  and  increases  roughly  like  for  higher  gain. 

This  nonlinear  behavior  for  threshold  is  expected  from  Eq.  (A-1)  where  ~  e^^  e'^. 
Figure  A-3  shows  ASE  pulses  for  p  =  10^  and  10^®,  g  =  1  and  various  p  values. 


The  experimentally  expected  ASE  pulse  should  be  calculated  by  integrating  I(t) 
from  Eq.  (A-1)  over  all  lengths  x  from  0  to  /  to  take  into  account  the  possibility  that 
spontaneously  emitted  photons  from  all  points  in  the  active  region  can  be  amplified  and 
contribute  to  the  measured  pulse.  Thus 


O 


(A-7) 


gives  the  expected  pulse  for  low  gain.  For  high  gain  on  the  other  hand,  a  simple 
spontaneously  emitted  photon  could  conceivably  depopulate  the  whole  system.  In  this  case 
different  runs  with  the  same  initial  conditions  would  produce  different  pulses  and 
Eq.  (A-1)  then  would  represent  the  result  of  many  experiments. 
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In  Figs.  A-4(a)  through  A-4(f)  we  show  a  surface  plot  of  |  aO  -x,t)|  as  a  function 

of  X  and  t  for  ^  =  0,5,10^0,40,80.  The  volume  under  the  surface  gives  |A(/,  t)  tmalf  and 
is  given  in  Table  A-3  for  each  value  of  p  along  with  the  functions  gi(ii)  =  and 

g2(H)  =  ■  for  comparison. 

The  firacdon  gi(p.)  describes  the  pulse  reduction  after  transmission  over  distance  1 
and  g20i)  describes  an  average  reductitm  after  transmission  over  different  distances  from  0 
to  L  This  represents  uniform  source  distribution  in  a  rod  and  pulse  emission  at  one  end. 

This  volume  (total  intensity)  calculation  does  not  follow  the  approximate 

dependence  obtained  in  the  last  column  on  the  right  in  Table  A-3  as  assumed  by  Baldwin 
(Ref.  A-3). 


Table  A-3.  Parameters  for  Describing  Pulses 


e-M^ 


1.3x10' 


Normaized 

Volume 


1.0 


5 

1.0x10® 

10 

3.0  X  10® 

20 

1.2  X  10 

40 

1.7x10"^ 

80 

2.0  X  10-2 

2.3  X  10“^ 


9.3  X 10"^ 


1.3x10-® 


1.5x10"® 


1. 


6.7  X  10~® 


4.5  X  10”® 


2.1  X  10"® 


4.2  X  10-^® 


1.8x10"®® 


1.0 


.2 


0.1 


0.05 


0.025 


0.0125 


1 


Figura  A->4<a*f).  A  two-dimensional  plot  of  I  A(x),  t)l  ^  the  pulse  Intensity  as 
a  function  of  time  for  different  values  of  the  attenuation  coefficient 


p  S  0,  2,  5.  10,  20,  40,  80,  Shown  In  (a),  (b),  (c),  (d),  (e),  and  (f). 


respectively. 
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